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ABSTRACT 

We show that when the gravitational force is correctly calculated in deal- 
ing with the vertical hydrostatic equilibrium of black hole accretion disks, the 
relationship that is valid for geometrically thin disks, i.e., c s /QkH = constant, 
where c s is the sound speed, fix is the Keplerian angular velocity, and H is the 
half-thickness of the disk, does not hold for slim disks. More importantly, by 
adopting the correct vertical gravitational force in studies of thermal equilibrium 
solutions, we find that there exists a maximally possible accretion rate for each 
radius in the outer region of optically thick accretion flows, so that only the inner 
region of these flows can possibly take the form of slim disks, and strong outflows 
from the outer region are required to reduce the accretion rate in order for slim 
disks to be realized. 

Subject headings: accretion, accretion disks — black hole physics — hydrody- 
namics 



Introduction 



The slim disk model is one of popular models for accretion flows around black holes 
and has been applied in recent years to quite many energetic astrophysical systems such as 
narrow line Seyfert 1 galaxies (e.g., Wang & Netzer 2003; Chen & Wang 2004), Galactic 
black hole candidates (e.g., Watarai et al. 2000), and ultraluminous X-ray sources (e.g., 
Watarai, Mizuno, & Mineshige 2001; Vierdaynti et al. 2006). 

Despite of its growing importance in the observational sense, we notice that some theo- 
retical problems regarding to this model in the fundamental sense have been ignored. In this 
paper we address one such problem, namely the inaccurate calculation of the gravitational 
force in dealing with the hydrostatic equilibrium in the vertical direction of slim disks. We 
work in the cylindrical coordinate system (r, z, (p) throughout. 
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2. Vertical gravitational force 

Obviously because the Shakura-Sunyaev disk (SSD) model (Shakura & Sunyaev 1973) is 
the first and still the most successful model of accretion disks, many simplifications made in 
this model are followed by the subsequent accretion disk models without rigorously verifying 
their applicability. The treatment of the disk's vertical structure is such a one. 

The very basic character of an SSD is that it is geometrically thin, i.e., everywhere in 
the disk the half thickness H(r) is much smaller than the corresponding cylindrical radius r, 
H(r)/r <C 1. This means that the averaged motion of disk matter in the vertical direction 
(if there is any) must be negligible comparing with that in the radial direction, and it is 
reasonable to assume that in the macro sense the disk matter is in the vertical hydrostatic 
equilibrium, i.e., in the vertical direction the gravitational force and the pressure force are 
balanced with each other, 

dp dip 

where p is the pressure, p is the mass density, ip is the gravitational potential, and the disk 
is already assumed to be steady and axisymmetric, d/dt — d/d(p — 0. In the SSD model the 
Newtonian potential ^(r, z) = —GMj\Jr 2 + z 2 was used, where M is the mass of central 
accreting object, and the one-zone approximation in the vertical direction was made; then 
in equation (1) dp/dz ~ —p/H, dip^/dz ~ GMH/r 3 = Vl^H , where Vt K = (GM/r 3 ) 1 / 2 is 
the Keplerian angular velocity, and a relation H = c s /Qk was obtained, where c s = (p/p) 1 ^ 2 
is the sound speed (e.g., Kato, Fukue, & Mineshige 1998, p. 80). 

The above procedure was somewhat improved by HOshi (1977), who approximated the 
Newtonian potential as ifa(r,z) ~ ifa(r,Q) + Q 2 <i z 2 /2 and assumed a polytropic relation in 
the vertical direction, p = Kp l+l l N , where K and N are constants, instead of the one- zone 
approximation. Then the disk becomes to have a vertical structure, and the vertical inte- 
gration of equation (1) gives H = \J2(N + l)c s /f2 K , where the sound speed c s = (po/po) 1 ^ 2 , 
with subscript representing quantities on the equatorial plane. 

Both the SSD model and HOshi (1977) reduced the differential equation (1) into a very 
simple but very useful relation: c s /Vty.H = constant. This relation was adopted in the 
slim disk model in the following way. First, the vertical hydrostatic equilibrium was still 
assumed, even though slim disks are not geometrically thin, i.e., they may have H < r or 
H ~ r. The argument for this was that these disks have quite large radial velocities, so that 
their configurations are quasi-spherical, i.e., the vertical motion of disk matter can still be 
safely neglected, and equation (1) holds. Second, although the pseudo-Newtonian potential 
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introduced by Paczyriski & Wiita (1980, hereafter the PW potential), i.e., 



# w (r, z) = -GMbh/(V^ + ^ - r g ) , (2) 

was widely used to simulate the general relativistic effect of central black hole, where M BH 
is the black hole mass and r g = 2GMbk/c 2 is the Schwarzschild radius, it was again treated 
in the form of Hoshi (1977), i.e., 

V> PW (r, z) ~ # w (r, 0) + Vi\z 2 /2 , (3) 

where the Keplerian angular velocity f2x = (GM/r) 1 ^ 2 /(r — r g ). Using equation (3) and 
keeping the assumption of polytropic relation, the relation c s /VL-^H = constant was refined 
by integrating equation (1), with the sound speed c s being defined either in terms of the 
pressure and mass density on the equatorial plane (Abramowicz et al. 1988; Wang & Zhou 
1999; Chen & Wang 2004), or in terms of the vertically integrated pressure and density, i.e., 
c s = (n/S) 1 / 2 , where IT = f" H pdz and £ = fH H pdz (Kato et al. 1998, p. 242; Watarai 
2006). 

To check the validity of the Hoshi form of potential, we draw in Figure 1 the gravitational 
force in the vertical direction dippw/dz in units of (c 2 /r g ) for varying values of (z/r) at a 
fixed radius r = 10r g , calculated from equations (2) (the solid line) and (3) (the dashed line), 
respectively. It is as expected that the Hoshi form of potential equation (3) is valid only for 
(z/r) < 0.2, while for (z/r) ~ 1 it greatly magnifies the vertical gravitational force comparing 
with the correct result according to the explicit form of the PW potential equation (2). 

In the following calculations we take the polytropic index iV = 3. Instead of the simple 
relationships (p/po) 1 ^ 3 = (p/po) 1 ^ = ^ — (z/H) 2 and8p /po = ^i^H 2 obtained using equation 
(3) (Kato et al. 1998, p. 241), the vertical integration of equation (1) using equation (2) 
gives 

JL] = (L\ = Vl+(^A) 2 -r g /r ~ y/l+(H/r)*-r s /r 

1 ~ r g/ r y/l+(H/r) 2 -r s /r 



PoJ \PoJ 



4^ = ^(1 - rjr) 
Po 



l-rjr 



+ (H/r) 2 - rjr 



(5) 



For a fixed r and a given value of H/r, we obtain values of (E/2p H) and (U/2p H) by 
vertically integrating equation (4), and value of (c s /VL K H) by using equation (5). We de- 
fine c 2 = n/E as in Kato et al. (1998), where detailed results were provided and can be 
quantitatively compared with ours. 
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Figure 2 shows T,/2p H (the dashed line), Il/2p H (the dotted line), and c s /VLkH (the 
solid line) as functions of H/r at r = 10r g . In the slim disk model these three quantities were 
all constant, for iV = 3 they were: £/2p # = 16/35, Tl/2p H = 128/315, and c s /Q K H =1/3 
(Kato et al. 1998, p. 242). It is seen from the figure that, however, these quantities are not 
constant, they take approximately their model-predicted values only for H/r <C 1. For slim 
disks with H/r ~ 1, the simple relation c s /VL-&H = constant is invalid. 



3. Thermal equilibria 

To reveal further the consequence of the correction of vertical gravitational force, we 
go on to study thermal equilibrium solutions of black hole accretion flows. We write the 
continuity, radial momentum, angular momentum, energy, and state equations of accretion 
flows in the form similar to that of Kato et al. (1998): 

M = -2nrY,v r = const. , (6) 



dlnr s dlnr "" E J_ H y\ + [ z / r f - rjrf^l + (z/rf 

M{Vtr 2 - jo) = 2nar 2 U , (8) 

Qtis = Qadv + Qr&d ' (9) 

n = n g „ + n», = ^ f (£_) (I.) iz + 1 T} f (l\\ z (10) 

/U^p J-h \PoJ \ToJ 3 J_ H \T J 

where M is the mass accretion rate, v r is the radial velocity, Q is the angular velocity, j is an 
integration constant representing the specific angular momentum (per unit mass) accreted 
by the black hole, and a is the viscosity parameter. Equation (7) is obtained by vertically 
integrating equation (8.2) of Kato et al. (1998) and specifying the PW potential equation (2). 
The viscous heating rate is Q^ is = MQ 2 fg/2ir, where f — 1 — j/Q^ 2 , j = jo/w, u = Q/Q-& 
is assumed to be a constant that is smaller than 1 (sub-Keplerian rotation) and is to be 
evaluated, and g = — dlnfi K /dlnr. The advective cooling rate is Q~ dv = £Mc 2 /2nr 2 (see Kato 
et al. 1998 for the detailed expression of £). The radiative cooling rate is Q rad = 32o-T 4 /3r, 
where T is the temperature on the equatorial plane, r = K, es poH is the vertical optical depth 
(e.g., Eq. [8.61] of Kato et al. 1998), and k cs ~ 0.34 cm 2 g _1 is the electron scattering opacity 
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that is assumed to be the dominant opacity source. In equation (10) T is the temperature, 
and T/Tq = (p/po) 1 ^ 3 as expressed by equation (4). For slim disks the radiation pressure 
dominates over the gas pressure, so the term Il gas in equation (10) can be dropped for the 
moment. 

To avoid unnecessary complicacy, we ignore the ram pressure term v r dv r /d\nr in equa- 
tion (7), and take dlntl/dlnr = —3/2 and £ = 3/2 from the self-similar solution (e.g., Wang 
& Zhou 1999; Watarai 2006). Then we have a set of six algebraic equations, i.e., equa- 
tions (5-10), which can be solved for six unknown quantities H, v r , c s , Q (or uj), po, and 
To, with given values of Mbh, M, a, j, and r. In our calculations we fix Mbh = 1OM , 
a = 0.1, and j = 1.83cr g (a reasonable value that is just a little less than the Keplerian 
angular momentum at the last stable orbit, 0, K r 2 \ 3rg = 1.837cr g ). 

Figure 3 shows thermal equilibrium solutions at a certain radius in the form of rh — H/r 
plane, where m is the accretion rate normalized by the Eddington accretion rate M^dd = 
QAttGMbh/cKcs. The thick lines represent the solutions obtained with the PW potential 
equation (2). For comparison, the thin lines draw the solutions obtained with the Hoshi 
form of potential equation (3), and accordingly equation (7) is reduced to equation (2.9) of 
Matsumoto et al. (1984). The dashed, solid, and dotted lines are for r = 5r g , r = 50r g , 
and r = 500r g , respectively. It is seen that for r = 5r g , the difference between the results 
obtained with two potentials is only quantitative. For r = 50r g and r = 500r g , however, 
the difference becomes qualitative. The most remarkable is that the thick solid and thick 
dotted lines have a maximum (not drawn in the figure), i.e., there exists a maximally possible 
accretion rate m max for slim disk solutions in the PW potential, while no such a m max exists 
for solutions in the Hoshi form of this potential (the thin solid and thin dotted lines). One 
might wonder whether H/r > 1, at which rh max appears, is physical. Our arguments are the 
following. First, H/r > 1 has indeed been obtained in many works on slim disks or other 
accretion disks (see, e.g., Figs. 5-10 of Peitz & Appl 1997; Figs. 1-4 of Popham & Gammie 
1998; Figs. 4 and 6 of Lu, Gu, & Yuan 1999; Fig. 5 of Chen & Wang 2004; Figs. 3 and 5 of 
Watarai 2006). Second, it does not matter for which value of H/r this m max really appears. 
What is important is that for each large radius there is an upper limit for m, beyond which no 
thermal equilibrium solutions can be constructed. By contrast, the previous understanding 
in the slim disk model (e.g., Abramowicz et al. 1988; Chen et al. 1995; Kato et al. 1998) 
was that any large m can correspond to a thermal equilibrium solution, as the thin solid and 
thin dotted lines in Figure 3 imply. This is because those authors either used the Hoshi form 
of potential equation (3) that magnifies the vertical gravitational force, or used the relation 
c s /QkH = constant that is resulted from equation (3). 

It is also noticeable that at small radii (r = 5r g in Fig. 3) there is no m max for thermal 
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equilibrium solutions. Figure 4 illustrates this qualitative difference between small and large 
radii. Figure 4(a) is for m = 10 and r = 5r g , Figure 4(b) is also for m — 10 but for r = 50r g 
{fn is larger than the corresponding rh max ). It is clear that in the former case the thermal 
equilibrium can be established, i.e., Q+ is = Q~ (= Q~ dv + <5 r " ad ), or Q+ is = Q^ v if Q'^ is 
totally negligible, can be realized for some value of H/r ~ 1; while in the latter case there 
are no thermal equilibrium solutions as Q^ is is always larger than Q~ . 

We derive an approximate analytic expression of the critical radius r crit , for r > r crit there 
is a m max and for r < r cr i t there is not. Equation (5) gives 4c 2 < ft^r 2 , and equation (7) 
gives (3/2)c 2 + VL 2 r 2 ~ ^k?" 2 , so we have (5/2)c 2 < Vt 2 r 2 . As seen from Figure 4, the 
condition that can ensure a thermal equilibrium solution at a radius is Q~ dv > <5v; s , i-e., 
(3/2)c 2 /r 2 > fi 2 /# ~ (3/2)fi 2 / since g ~ 3/2. Therefore, the criterion for a thermal 
equilibrium solution to exist is expressed as 

f = 1 ^— < - . (11) 

For j = 1.83cr g in our calculations, the critical value / crit = 2/5 gives r crit = 16.4r g , which 
agrees well with the numerical value r crit = 18.8r g in Figures 3 and 5. The reason why / is 
of crucial importance is that since Q^ is is proportional to /, for small radii / is small, and 
Qvis can be balanced by Q~ dv for any to; while for large radii / is large, there must be an 
upper limit for to, beyond which Q^ is would be too large to be balanced by any cooling. 

To see where the position of slim disks is according to our understandings, we present 
in Figure 5 a united description of thermal equilibrium solutions of optically thick accretion 
flows around black holes in the form of m—r plane. The solid line in the figure draws m max for 
each radius, above which no thermal equilibrium solutions exist at all. The rest of the plane 
is further divided into three regions based on the local stability analysis. We adopt a simple 
thermal instability criterion: (dQ+ is / dT)% — (dQ~ / dT)^ > 0, which is practically valid for 
moderately large-scale perturbations (Kato et al. 1998, p. 306). The criterion results in: 
5 = 2-5(3- 6/ adv + 8/3/ adv > 0, where (3 = IT gas /IT, and / adv = Q; dv /Q+ S - The dotted and 
dashed lines both draw 5 = 0, but correspond to f3 — 2/5 and / adv = 0, and to / adv = 1/3 
and (3 — 0, respectively. The region below the dotted line has 5 < and is for stable SSDs, 
which are radiative cooling-dominated (/ adv — > 0) and gas pressure-supported {(3 > 2/5). 
The region between the dotted and the dashed lines has 5 > and is for unstable SSDs, 
which are radiation pressure-supported ([3 < 2/5) but not yet advective cooling-dominated 
(/adv < 1/3). These two regions were already known in the literature. The region between 
the dashed and the solid line has 5 < and is obviously for slim disks, which are advective 
cooling-dominated (/ adv > 1/3) and radiation pressure-supported {(3 — > 0), and are stable. 
What is new is that, because of the limitation of rh max , accretion flows can possibly take the 
form of slim disks only in the inner region, i.e., in the region r < r crit . 
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Figure 5 may have interesting implications. If the accretion rate m of an accretion flow 
is sufficiently small at large radii (i.e., below the dotted line), then the flow can behave as 
an SSD throughout (here we don't want to discuss the inability of SSD model for the inner 
region of a black hole accretion flow, such as for the flow's transonic motion). If, however, 
rh at large radii is in the unstable region or in the no-solution region, then it seems that the 
only possible way for accretion to proceed is that the flow loses its matter continuously in 
the form of outflows. Such outflows must be so strong that m keeps being below the dotted 
line all the way till the inner region r < r crit . In the inner region there are two possibilities 
for the accretion flow: either m is still below the dotted line and the flow keeps being an 
SSD, or m is in the unstable region. In the latter case the flow may undergo a limit cycle, 
i.e., oscillating between the SSD state and the slim disk state in the vertical axis direction 
of Figure 5. Such a limit cycle behavior has been investigated extensively in the literature 
(e.g., Szuszkiewicz & Miller 2001 and references therein), and is the likely way for slim disks 
to be realized. A related remark is that provided m in the outer region of accretion flows 
is not sufficiently small, outflows seem to be unavoidable, as already observed in many high 
energy astrophysical systems that are believed to be powered by black hole accretion. 



4. Discussion 

We have shown that when the gravitational force in the vertical direction of black 
hole accretion disks is correctly calculated, the relationship c s /QkH = constant, which is 
valid only for geometrically thin disks, does not hold for slim disks; and that there exists 
a maximally possible mass accretion rate for each radius in the outer region of optically 
thick accretion flows, so that only the flow's inner region can possibly take the form of slim 
disks, and outflows from the outer region must be produced in order for slim disks to be 
realized. We stress that only one change has been made in obtaining these results, i.e., 
using the explicit form of the PW potential, instead of the Hoshi form of this potential, to 
calculate the vertical gravitational force, while all the assumptions, equations, and methods 
for solutions are kept exactly the same as in the slim disk model (e.g., in the excellent book 
of Kato et al. 1998). 

All our results here are based on a local analysis. Though the similar local analysis was 
often used in the literature (e.g., Abramowicz et al. 1995; Chen et al. 1995; Kato et al. 
1998), it is worth checking our results by constructing global solutions of original differential 
equations for black hole accretion flows, like what was done by, e.g., Chen & Wang (2004), 
but with a revised vertical gravitational force. We plan to do this in a subsequent work. 

A natural question is why we concentrate our attention here only on slim disks and 
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do not touch another type of black hole accretion flows, namely optically thin advection- 
dominated accretion flows (ADAFs, Narayan & Yi 1994; Abramowicz et al. 1995), which are 
also geometrically not thin. We think that the ADAF model is also likely to suffer the same 
problem in the vertical direction as the slim disk model does, because the problem is resulted 
from a purely hydrodynamic consideration and is related only to the geometrical thickness of 
the flow. In particular, the usage of relation c s /VL-^H = constant for ADAFs is questionable. 
What is different is that for optically thin flows the radiation processes are more complicated 
and the vertical integration is not as easy to perform as for slim disks. In addition, ADAFs 
are known to correspond to very low accretion rates and to have a maximally possible 
accretion rate at each radius (e.g., Abramowicz et al. 1995), so the problem may have no 
impact on ADAFs in respect of the accretion rate. However, we still wish to make one more 
comment. Many 2-dimensional or 3-dimensional numerical simulations of viscous radiatively 
inefficient accretion flows revealed the existence of convection-dominated accretion flows 
instead of ADAFs (e.g., Stone, Pringle, & Begelman 1999; Igumenshchev & Abramowicz 
2000; McKinney & Gammie 2002; Igumenshchev, Narayan, & Abramowicz 2003). According 
to our results here, the 1-dimensional ADAF model might have hidden inconsistencies in the 
vertical direction, so that ADAFs could not be obtained in those multidimensional numerical 
simulations. 

We thank Li Xue for beneficial discussion and the referee for prompt and helpful com- 
ments. This work was supported by the National Science Foundation of China under Grants 
No. 10503003 and 10673009. 
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Fig. I. — Vertical gravitational force dijj/dz for varying z/r at r = 10r g , calculated using 
the explicit form of the PW potential (the solid line) and the Hoshi form of this potential 
(the dashed line), respectively. 
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Fig. 2. — Variations of T,/2p H (the dashed line), U/2p H (the dotted line), and c s /VL-^H 
(the solid line) with H/r at r = 10r g . 
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Fig. 3. — Thermal equilibrium solutions with the Hoshi form of the PW potential (thin 
lines) and with the explicit PW potential (thick lines) at r = 5r g (dashed lines), r = 50r g 
(solid lines), and r = 500r g (dotted lines). 
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Fig. 4. — Variations of Q^ is (the solid line), Q~ dv (the dashed line), Q~ d (the dotted line), 
and Q~ dv + Q-rad (the dot-dashed line) with H/r. (a) is for m — 10 and r = 5r g , (b) is for 
m — 10 and r = 50r g . Numbers attaching the vertical axis are in units of erg cm _2 s _1 . 
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Fig. 5. — Distribution of thermal equilibrium solutions. The solid line is for m max , the 
dashed line for / adv = 1/3, and the dotted line for (3 = 2/5. 



